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Abstract 

We present a calculation of the width for Xcj —>■ e^e“ decay. The amplitude of the process is 
computed within the NRQCD framework. The leading-order contribution is described by two terms 
associated with the two different integration domains in the electromagnetic loop describing two- 
photon annihilation of the heavy quark-antiquark pair. The corresponding operators are defined in 
the framework of NRQCD. The matrix element of one of these operators describes a configuration 
with an ultrasoft photon and can be associated with the higher Fock state of the heavy meson. In 
order to compute this contribution we use the heavy hadron chiral perturbation theory. We obtain 
that this contribution is numerically dominant. The obtained estimates for the decay widths of the 
Xci and Xc 2 states are 0.09 eV and 0.07 eV, respectively. 
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1 Introduction 


The leptonic decays of C-even charmonium states into a lepton pair have a very small branching ratio 
because they can only occur via a two-photon intermediate state XcJ 7 * 7 * —>■ e+e”. However with the 
high-luminosity BEPC-II e“'"e“-accelerator operating on the charmonium energy region such measure¬ 
ments of direct production cross section e“'"e“ —>■ XcJ become feasible. The study of the mechanism of 
the production of C-even quarkonium states is especially interesting in view of the production of higher 
resonances such as the exotic charmonium-like state X(3872). 

The decays XcJ —t l~^l~ have already been studied long time ago in Ref. [^. The authors consid¬ 
ered the quarkonium description and a vector dominance model (VDM) in order to describe the decay 
amplitudes of charmonium states with J = 1, 2. It was found that a naive quarkonium description is prob¬ 
lematic because of infrared logarithmic divergencies arising in the integrals describing the quark-photon 
annihilation loop. Such divergencies indicate that the corresponding contribution is also sensitive to long 
distance physics. The corresponding integrals in Ref. have been regularized by introducing the binding 
energy — 2mQ ~ 500 MeV. The numerical estimates obtained in this way give very small value of the 
widths, smaller thn bounds derived from analyticity and unitarity in the same work. Probably, a more re¬ 
alistic estimate was obtained using generalized VDM which yields a larger numerical value for the widths, 
consistent with the unitarity constrains: r[xci —t e+e“] ~ 0.46 eV and r[xc 2 —t e+e“] ~ 0.014 eV [^. 

Recently the decay rate of the Xci state was again estimated in using the VDM approach, with 
result r[xci —t e+e“] ~ 0.1 eV. The short distance contributions describing a configuration with highly 
virtual photons were considered in this framework as unknown contact vertices giving rise to a theoretical 
uncertainty. 

The aim of the present work is to provide a more systematic description of the decay amplitudes for 
XcJ 7*7* —t e+e“ process using the NRQCD factorization framework |^[^, see also review and 
references theirin. This technique allows one to perform a systematic description of heavy quarkonium 
states using the small relative velocity of heavy quarks and the small QCD running coupling at short 
distances. Within this framework we associate the IR-divergencies found in with a specific quark- 
photon operator which describes a configuration with an ultrasoft photon. The matrix element of this 
operator describes the overlapping with the higher Fock state which consists of heavy quark-antiquark and 
photon. This allows us to perform a systematic separation and description of the different contributions 
relevant in a leading-order expansion in small velocity v. 

Our paper is organized as follows. In Sec. [2 we briefly describe our notation and provide definitions 
of various quantities used in the following. Sec. ^is devoted to the investigation of the one-loop integral 
which describes the leading-order contribution. In this section we establish the dominant regions and 
provide a description of the amplitude within the NRQCD factorization framework. Sec. is devoted 
to the calculation of the ultrasoft photon matrix element in the heavy hadron chiral effective theory 
(HHxPT). Furthermore, a estimate of the decay rates is given. In Sec. we briefly summarize our 
results. 


2 Kinematics and notation 

Let us start from the description of the decay kinematics Xcj{P) e“''(Zi)e“(/ 2 )- The initial state 
momentum can be written as 

P = M^uj, uj^ = I, (I) 

where w denotes the charmonium velocity. In the following, we consider the charmonium rest frame 
which implies 

w = (I,0). (2) 

The small relative velocity of heavy quarks in the bound state is denoted as v. Neglecting lepton masses, 
the lepton momenta can be written as 
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h=M^-, h=M^-, (3) 
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where n and n denote the light-like vectors which satisfy {n • n) = 2. Any 4-vector can be expanded 

as 

V>^ = {V-n)^ + iV-n)'^ + V^, (4) 

where V± denotes the components transverse to the light-like vectors : (Vj_ • n) = (Vj_ • h) = 0. Similarly, 
one can also write a decomposition 

= iV + V-^, (5) 


where Vy denotes the component which is orthogonal to the velocity uj: (w • VV) 
we assume that in the rest frame 


n n 



0. In the following 


( 6 ) 


The momenta of the heavy quark and antiquark with mass m which form a quarkonium state can be 
written as 

+ P2 = ^P-A, (7) 

where the relative momentum A satisfies 


(A-w) = 0, A2 = -A2. 


( 8 ) 


The heavy quarks which create a bound state are non-relativistic, implying that the relative velocity 
V A/m is quite small: z) ^ 1. 

The power counting rules of NRQCD has been established in [^|^. Following these arguments we 
assume that the mass m is large enough and that the most important scales such as mass m, typical 
three-momentum of the heavy quark ^ mv and its typical kinetic energy ~ satisfy 

{mv^Y (mv)^ ^ m?. (9) 

Integrating out the modes with hard momenta ph m^ one passes onto the effective theory NRQCD 
which describes the modes with the soft momenta Ps ~ mv. If the scale mv Aqcd then one can 
integrate over the soft region together with potential gluons with momenta |^[7- 

Po ~ mv'^, p ^ mv, (10) 
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After this one obtains a new effective theory which is known as potential NRQCD (pNRQCD). For a 
more detailed information about these effective theories see Ref. and references therein. 

The charm quark mass me ~ 1.5 GeV is not large enough compared to Aqcd therefore in this case one 
can only factorize the effects at momentum scales of order me- However, in the QED sector one can also 
consider the possibility to integrate over the soft region too. As we will show further on, such situation 
is relevant for the factorization of the electromagnetic loop describing the transition cc —>■ 7 * 7 * —>■ e“'"e“. 

After factorization of the hard contribution, the nonpertubative QCD dynamics is described by the 
matrix elements of appropriate operators defined in an effective theory. In the following, we will need 
the following set of NRQCD operators which describe the matrix elements between the charmonium and 
the vacuum states: 

( 11 ) 

Oi^Po) = ( 12 ) 




0^^{^P2)=xl 
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7t:7t 
/ 3 ) ; 




(13) 

(14) 


where the covariant derivative -|- gA^^, D -y = D-y — D-y- Furthermore, 

four-component fields Xio to describe the soft quark and antiquarks within the 
These fields satisfy 




we use the covariant 
NRQCD framework. 

( 15 ) 
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Using Eq.([^, one can show that the operators in Eqs. can be reduced to the set of well-known 

non-relativistic operators constructed from two-component Pauli spinors. 

The matrix elements of these operators are well known in the literature and can be written as 


(0|V'i7TX-l>/M=4 (16) 

(O'CS^)), (17) 

{0\0{^Po)\xco)=i{OCPo)), (18) 

(0|(!l"(3Pi)IXci)=*e^ (OCPo)), (19) 

( 0 | 0^^CP2) | xc 2 ) = ief {O^Po)) . ( 20 ) 


The constants (©(^S'l)) and (©(^Pq)) are related to the value of the charmonium wave functions at the 
origin 


(0(35i)) = ^ i?io(0), 

(0'(3^i)) = P2 o(0), 


( 21 ) 

( 22 ) 


(OCPo)) = v^^/ 2 iw^ 



(23) 


where Rnijr) is the radial part of the Schrodinger wave function and R'ni{f) denotes its derivative. The 
rhs of Eqs.(18l-20) depends on of the same constant (Oi^Po)) due to the spin symmetry of the leading 
non-relativistic action li). The polarization vectors eVA), eJ(A) and e“^(A) correspond to spin-1 and 


spin-2 charmonium states, respectively. They are normalized to satisfy 


pa pp 

^ e5c(A) {e^(A)}* = (24) 

A ^ 

with X = {J/ip, Xci} and 

e“^(A) |e“ ^ (A)| = -|- -Map'Mpa' — , (25) 

with Map = -gap + PaPp/M^^^. 

The decay amplitudes Xcj —t e+e” are defined as 

(e+e“; out| *n; Xcj) = i{2TT)^5{li +I 2 - P) Aj, (26) 


with 

A,J = UnTjVn Tj, (27) 

and where and Vn denotes the spinors of the massles lepton and antilepton, respectively 


u„ 




(28) 


and 

hi = r2 = (29) 

The leading-order contribution to these amplitudes arises from the annihilation of heavy quarks into 
two photons which create the outgoing lepton pair, see Eig. If the one-loop integral is dominated by 
the hard region where both photons and heavy quark are highly off-shell then one can expect that such 
process can be described within the NRQCD approach and the amplitude can be factorized into hard 
and soft parts. In the next section we consider this possibility in more detail. 
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Figure 1: One-loop diagrams describing the annihilation of Xcj into an e“'"e pair. 


3 Factorization of decay amplitudes in NRQCD 


The leading-order in Og diagrams describing the e“'"e“ decay of C-even charmonia are shown in Fig. 
These one-loop diagrams are constructed from the photon, lepton and heavy quark (double lines). The 
diagrams in Fig. [^can be computed in the heavy quark mass limit m —> oo, performing an expansion 
in the small parameter A/to ~ v. Let us to start from a naive guess that the dominant contribution is 
only provided by the hard region where the loop momentum ~ to, and therefore all propagators are 
far of off-shell. The leading-order contribution in 1 /to is provided by projections onto the leading-order 
operators Oi^Pj) described in Eqs.(18)-(20|. The technical details are well known in the literature, see 
e.g. [^. The resulting expressions can be presented as 

Unjaih - 


Ai=e';^i{0{^Po))e^V2 / dk j 


A 2 — Sy 


i{Oi^Po))e^ / dk 


'(M 

1 

1 _ 1 

[P] 

1 

to 

'^n^a. (/l 


1 - 1 

1 

to 

■ [P] 

1-1 

1 
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^Tr [V^r°‘^f^{k)] 


(30) 


(31) 


where the square brackets for the propagators denote the standard Feynman prescription [A] ^ = 
[A -|- ie] . The corresponding contribution to the amplitude XcO e+e“ vanishes and therefore is 
suppressed by a power of v and will not be considered it in this work. 


The total structure of the integrands in expressions (30) and (31) can be divided into the lepton and 


heavy quark parts. The lepton part has Un ■ ■ - Vn in the numerator and includes the photon and lepton 
propagators in the denominators. The heavy quark part is given by Tr[7^ir“^^(fc)]. We introduced the 
projections Vj onto charmonium states 


' Pi '' = + ‘/^)(7t7t - 7t7t)75, 


’^2 — (1 + </’)7t- 

The expression for r“^'^(fc) reads 

^ {pDQik) + DQ{k)p] + 

with 

^2 

{7^(to</> — ^ + 771)7“ -I- 7“(|^ — TOt/ -I- to)7^} 


Dq = 




iiieeqf 


— 2m{kuj) — 

— 2m{kuj) — A^ 


(32) 

(33) 

(34) 

(35) 


-I- ^ m ) x °‘ — 7 “(^ — TO(/’ -I- to) 7 ^} , ( 36 ) 

k"^ — 2m{kui) — A^ 
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where eg is the charge of the heavy quark (cc = 2/3). The small squared relative momentum ^ (mv)'^ 
which appears in the heavy quark propagator provides an IR-regularization and can be neglected if it 
is not required. With this regularization the traces and loop integrals are computed in four dimensions 
with dk = d^k/{2'nY. 


The expressions (301 and (31) have been obtained by expanding the heavy quark fields in position 
space 

l + yd + ^ilpT 
2m 


c{y) ~ 


V’^(O) 


(37) 


and projecting the soft quark fields xtu and onto leading-order operators (121-( 14). The terms ^ y 


(arising from the multipole expansion of the soft quark field arguments ) lead to the expansion of the 
integrand with respect to small relative momentum A giving the contribution D'q. The terms proportional 
to ~ 2 mgive the contribution with Dq. The evaluation of the integrals in Eqs. (30) and (31) gives 


A 2 = UnT2V„ i{0{^Po)) -^e.Q ^21nA^/m^ -b ^ (1: 


m 

(38) 

(In 2 — 1 -b itt) j . 

(39) 


These expressions are in agreement with the results obtained in Ref. . We obtain that both amplitudes 
depend on the large logarithm ~ InA^/m^ which is sensitive to the soft scale A^. This shows that the 
starting assumption about one dominant region fc ^ m is incorrect. There must be at least one more 
domain where some propagators in the loop integral are soft. One can expect that the additional region 
is associated with the configuration when one of the photons is soft. In this case the propagator of the 
heavy quark is also soft and the hard configuration is described by the tree level subdiagram describing 
the annihilation cc —>■ e“'"e“ through one photon. 

In order to get an idea about the explicit definition of this region it is useful to investigate the integrals 
of diagrams in Fig. [^within the threshold expansion technique worked out in Ref. |^. According to this 
analysis the threshold kinematics is described by the following regions 


hard : kn 


m, 


soft : k,, ^ mv. 


potential : ko ^ mv^, k ~ mv, 
usoft : ku mv^. 


(40) 

(41) 

(42) 

(43) 


The same regions can also be considered for the photon with momentum P — k. These regions can be 
associated with the fields appearing in the effective Lagrangians, see e.g. Ref. [^. 

According to the threshold expansion prescription an integrand is expanded in each domain to a 
required accuracy and the resulting integral is computed in dimensional regularization. A detailed analysis 


which enters in Eq. (30) 


of the full expressions in Eqs. (30) and (31) is quite similar. To be definite let us consider the integral 


J = dk 


Unlaill - k)-ipVf, 


[k-hf [P] (k-pf 


-Tt [P 






(*)]• 


(44) 


Keeping the denominators of the heavy quark propagators in unexpanded 

{pi — k)"^ — m? =^ip-|-A — fc'j — w? = — Po^o + 2(fc • A) — A^ -|- ^Pq — (45) 


(p2 -fc)^-m^ = ( ip-A-fc') = fc2-Pofco-2(fc-A)-A^-bip^-m^, (46) 
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where Pq ''tt, Pq/4 — ^ (mw)^ , A ^ mv. In the hard region, the small scalar products with A 

and the term Pg /4 — ^ m? can be neglected resulting in 


P ± A — fc 1 — — (kP), 

- h 


(47) 


which appear in the expressions (35) and (36) (up to small regularization term A^ ). From dimensional 
counting one immediately finds 

(48) 

where F denotes the Dirac structure. Computing the hard integral Jh in dimensional regularization one 
finds the IR poles 1/e. These singularities must cancel in the sum with other contribution. 

Expanding the integrand (44) in the soft region yields 


Js — dk 


UnlahlpVn 

—2 {kli)\[k'^][ArrP] 4 


Tr [Pi^.r:'''^(fc)] , 


where jg given by (34) with 

[DqI - 

1 1 


—2{kui)] 


- 1)7“ + 7“(1 - (/)7^} 


(49) 


(50) 


n'A* 


{7V7“+7“7V} + 


1 




S 2m [-(A:a;)] ' ' ' ' ' J ' 2m [-(fcw)]^ 

Calculating the trace and performing the contractions in the numerator results in 

1 MnElUft 


{7^((/ + 1)7“-7“(1-</)7^}. 


Js ~ — -UnTiVn I dk - „ 

m3 J [fc2] [-(kuj)f 


m'^ 


(51) 


As the integral in (51) is scaleless it therefore vanishes in the dimensional regularization, i.e. Jg = 0. 


In the potential region (42), the expansion of the heavy quark propagator reads 

2 


^P ± A — kj —w? ~ Pg /4 — m^ — PoA:g — 

p 


{k±A). 


The computation of the corresponding integral then yields 


Jp ~ —UnTiVfi [ dk 


1 




Po/4 - - Pg/cg - + A^ 


(a a -a) 




m"’ 


(52) 


(53) 


However the poles in ko in the integrand of Eq.(53) lie in the same imaginary half-plane and therefore the 
integral over fcg vanishes. This observation is also true for the higher order contributions in v appearing 
from this domain. We can therefore conclude that the potential region cannot contribute in this case. 

In the ultasoft domain (431, the heavy quark propagators are expanded as 

2 


fc — -P ± A ] — m^ 


Performing the expansion of the integrand one gets 


~ PgV4 - m2 - Pgfcg - A^. 


JII, 


— UnTiVn I dk 


[F] Po74 - m2 - Pgfcg - A2 


.1 . 


(54) 


(55) 
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This integral has the same scaling behavior ^ m. as the hard integral Jh in Eq.(48). One can also see 


that the integral in Eq.(55) is UV divergent. The similar analysis can also be carried out for the second 


photon with momentum k — P. Therefore we conclude that the exact integral must be given by sum 


J — Jh T Ju 


(56) 


where Jus denotes the contributions from the both ultrasoft domains. This conclusion can be checked by 
explicit calculations. A similar conclusion for the two-photon diagrams in Figjl] has also been obtained 
in Ref. 


11 


Guided by this consideration we suggest that the additional relevant domain is described by the 
ultrasoft region. In order to find the description of the appropriate operator in the effective theory one 
has to integrate out hard and soft photons and leptons. After that the description of QED sector includes 
only collinear leptons and ultrasoft photons. The integration of the soft photon with the lepton and quark 
must be described in the framework of the effective theory. 

Within the above picture the factorization of the decay amplitudes can be described as a sum of two 
contributions 


Aj = UnTjVu cA) i{oCPo)) 

+ {e+e-\UQ)Yl{GhlYnmh{0) \xcj) • (57) 

The first term on rhs of this equation corresponds to the hard domain with the hard photons, 
denotes the corresponding hard coefficient function. 

The second term on rhs of Eq. (571 corresponds to the domain with the ultrasoft photon. The operator 
C>'^(^S'i) is defined in Eq.(ll). The outgoing collinear leptons are described by fields and which 
defined as 


^„(a;) = i’cix)- 


in ^ '4^c{x)- 


(58) 


The photon Wilson lines and Yn describe the interaction of the ultrasoft longitudinal photons with 
the energetic lepton and antilepton and read 


Y^ (0) = Pexp 


{“i 


ds n ■ B^^{sn) > , 1^(0) = Pexp < — 




ds n ■ B'^^{sn) 


(59) 


where denotes the ultrasoft photon field. The appearance of these Wilson lines is related with the 
fact that in a general gauge the tree level diagram with attachments of n ■ 5“'* photon to the collinear 
field describing the outgoing leptoij^ are resummed to the P-ordered exponents 


tpc 




1 + en-B^ 


1 


i{n ■ D) 




The leading-order hard coefficient function Cj is defined by the diagram in Figj^and reads 


^7 — 2 ^Q' 


(60) 


(61) 


The soft and collinear modes in the effective action describing the QED sector are decoupled. This 
property is well known in the soft-collinear effective theory, see e.g. Refs. 12 14 . This allows us to 

(62) 


contract the lepton fields in the second matrix element in Eq. (571 and rewrite it as 


with 


(e+e-| en(0)r„n0)7l>^f.(0)ef.(0) \xcj) = Unll^h 

OA^Si) = Y^(0)Yf,(0)O^(^Si). 


a^(^Si)ixcj}, 


(63) 


The presence of the soft scale in Eqs. (38) and (39) can be explained by the contribution with 
ultrasoft photon. Therefore in order to find the hard coefficient functions we have to perform the 

^We assume electrical charge is measured in proton units (positron is particle and electron is antiparticle) that allows 
to use the same notation for the covariant derivative and Wilson lines as in QCD. 
















Figure 2: The hard one-photon exchange diagram. 


matching onto the configuration described by Eq.(57). For that purpose we need to compute the ultrasoft 
matrix element (62) in the effective theory. 

The interaction of ultrasoft photons with quarks are described within the pNRQED. The ultrasoft 
photons have momentum p ^ mv^ so that photon field scales as 


TJUS 


mv^. 


The scaling of the quark fields reads 

V'a; -- ditjju. 


{mv)ip^, doipcj ~ 


(64) 


(65) 


Using this counting one finds 

fj(j) o%^Pj) ~ m-^{mvf. ( 66 ) 

At the same time 

~ ~ m-‘^{mvf. (67) 

However the pure quark operator 0°’(^S'i) is C-odd and therefore it cannot contribute to the matrix 
element with a C-even charmonium state 


O"('5i)IXcj)=0. 


( 68 ) 


In order to obtain a nontrivial contribution one needs to consider at least one interaction of an ultrasoft 
photon with the quark in pNRQED. We only need the two-particle sector describing the electromagnetic 
interactions of quarks (in rest frame uj = (1,0)) 


= J d^x (^iuj ■ d + ^ 


= J d^x (^)^o 


X ■ dx eeq oj ■ B^^{xo) H- eeQB^^{xo) ■ idx 


’ipu^ix), 


(69) 


(70) 


and analogous contributions with antiquark fields. The arguments of the ultrasoft photon field are 
expanded because the space components of the quark fields varies at x ~ Ijmv, the measure scales 


as dxo ~ 1/r 




{mv) 


-3 


With these rules one finds that 


and 


leading-order term (691 provides the soft quark propagator 

A<^(fc) = 


(wfc) — k"^/2m + ie 

A nontrivial contribution to the matrix element (0|... \xcj) can be obtained from T-product 


The 


(71) 


(72) 


which is of the same order as the hard contribution in Eq.(66|. Calculation of this T-product gives 
diagrams shown in Eig|^ The dashed lines can be associated with the collinear leptons or equivalently 
with the ultrasoft Wilson lines (591. These diagrams induce a mixing of the operators 01((^S'i) and 
0{^Pj) due to electromagnetic interaction in the framework of pNRQED. 

In order to perform the matching onto operators according to formula (57) one has also to compute 
the contribution of the diagrams in Figj^ The simplest way to proceed is to follow the same technique 
as we used above for diagrams in Fig{^ 
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Figure 3: The diagrams generated by the T-product (72) in pNRQED. The crossed circle denotes the 
vertex of the operator the dashed lines represent the Wilson lines associated with light-like 

directions n and h. The small crosses on the dashed lines show all possible attachments of the photon. 


Let us consider Xc 2 as initial state. In this calculation we set Pq = 2m and only keep the relative 
momentum A. Then the sum of all four diagrams gives 


2 

{e+e-\C^T{0;CSi),Cr[Bn} \Xc2) = i {0{^P^)) J. 


2m 


(73) 


Computing these diagrams we project the soft quarks fields on the operator Oi^P^) and substitute the 
corresponding matrix element which gives the factor e’^^i the coefficient 4 arises from the sum 

of the four diagrams shown in Fig|^ the ultrasoft loop integral reads 


Jus = {-i) j dk 


^ ^—{ojk) — A2/2m 


(74) 


This integral coincides with the ultrasoft integral of Eq.(|55|) obtained within the threshold expansion 
approach up to term Pq/4 — m^ which vanishes because we set Pq = 2m. The integral in Eq.(74| is 
UV-divergent and we use dimension regularization D = 4 — 2e in order to compute it. The result reads 


Jus - 


tD/2 


(271-) 


D 


A" 

m^p 


-2e 


(75) 


where ^p is the factorization scale. The 1/e pole is the UV-pole which describes UV-mixing of the 
operators 0;)(^S'i) and O'^i^Pj), schematically 


[ o ^( 35 i )]^ = + Zj OCPj), 


(76) 


where [0]p on the Ihs of Eq.(76) denotes the renormalized operator. Eurthermore, Zj ~ e^/e is the 
corresponding renormalization constant. Assuming MS-subtraction scheme one finds 


1 A^ 

['Ais] R — e n ^71 . 

^ 47r2 m^ip 


Hence we obtain 


{e+e-\C^T{ 0 ;CSi),Cr[B^^]}\Xc 2 )ii = UnT 2 Vni{OCPo)) 2In 


A^ 

mfip 


(77) 


(78) 


The soft matrix element for the Xci can be computed in the same way, resulting in 


{e+e-\C^T{O^^CS,),Cr[B^^]}\Xci)n = UnTiVn i{OCPo)) 2^2In 




mfip 


The hard coefficients are given by 


ciJ) - ~ Tills (0| O^i^Si) \xcj) 

- Ur^TjVn (0| 0-{^Pj) \Xcj) ' 


(79) 


(80) 
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Figure 4: An example of the diagram in pNRQCD describing the ultrasoft matrix element of Eq. (62|. 
The doted lines denote potential gluons with momenta given by Eq.(42). 


where the expressions for Aj are given by Eqs.(38) and (39). The important check of the factorization 


formula (57) is the cancellation of the ultrasoft scale in the expressions for obtained from Eq. (801. 


Substituting the computed expressions in Eq.(80) we obtain 


= v/21n^ 


2 ’ 
F 


(7(2) _ ^ 2 
<-77 - ^ 3^0 


In ^ + I (In 2 - 1 
3 


■ ztt) 


(81) 


(82) 


These expressions are the main result of this section. We observe that the soft scale cancel in Eqs.(81) 


and (82) as it must be. Hence the factorization formula described by Eq.(57) describes properly the 


ultrasoft region of the one-loop diagram. 

( 2 ) 

The coefficient function has an imaginary part which originates from the two-photon cut. Such 
mechanism can not work for Xci state, therefore is real. 

The hard coefficient functions depend on the factorization scale pip. Therefore 


Aj = UrAjVn (0| \Xcj) (mf), 

and the independence of the amplitude Aj on fip yields the evolution equation 


^riFj'Cn^ {(^( .^o)} MF ^77^ (/^f) — '^n MF 


\o;c^s,)\xcj) 


(83) 


(84) 


The solution of this equation depends on the initial condition defined at some scale fiQ. Performing 
numerical estimates one has to fix a value of this scale. By derivation this scale separates the hard region 
(two hard photons) from the ultrasoft region (hard and ultrasoft photons). Therefore it is natural to 
associate this scale with the virtuality of the ultrasoft photon and to set /tq to be of order 300 — 500 


MeV. Then the matrix element of the operator OA^Si) on the rhs of Eq.(83) describes only the ultrasoft 


nonperturbative contribution which can be only estimated within some iow-energy effective theory or 
model. Similar to the well known color octet mechanism, the operator can also be associated 

with the electromagnetic mechanism. The corresponding matrix element can be interpreted as an overlap 
with the higher Fock state \QQx) which includes a dynamical photon while the matrix elements of 
the operators Oi^Pj) describe the coupling to the dominant quark-antiquark state. Therefore the full 
description of the leptonic decay requires a knowledge on the subleading structure of the quarkonium 
state. 

In the large mass limit m —>■ oo one can consider a specific situation known as the Coulomb limit when 
the binding energy is larger then the typical hadronic scale E ~ mv^ Aqcd- In this case the strong 
coupling is quite small aa(mv) ^ v and ultrasoft contribution can be estimated within the pNRQCD. 
Then one has to compute the diagram as in Figj^resumming the interactions with Coulomb gluons. Such 
calculation has been carried out for the radiation function in Ref. 15 . Perhaps, such calculation might 


also be interesting here in order to get an idea about the relative value of this matrix element in the 
Coulomb limit. In present paper we will obtain an estimate of the ultrasoft matrix element using the 
so-called heavy hadron chiral perturbation theory (HHyPT) framework in the next section. 
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4 Phenomenology 


4.1 Calculation of the ultrasoft matrix element in the heavy hadron chiral 
perturbation theory 

In order to provide a numerical estimate of the decay rate we need to estimate the ultrasoft matrix 
elements 

(0| 0^(35i) \xcj) (85) 

which describes an overlap with the higher Fock component of the charmonium state XcJ in which a 
dynamical photon is present. One can expect that the soft photon has already quite large wavelength 
and therefore it interacts with heavy charmonium as with a point-like source. Then it is natural to 
expect that the relevant dynamical degrees of freedom in this case are associated with mesonic fields 
and the corresponding low energy dynamics is described by the most generic effective action compatible 
with the symmetries of NRQCD. Such an approach is known as heavy hadron chiral perturbation theory 
in Refs. 16 17 for the heavy-light mesons and then generalized on quarkonia in Refs. 18 -20 . This 


framework can also be used for the calculation of the matrix element in Eq.(85). 


For our purpose we need only the electromagnetic sector of the HHyPT described by the effective 
action which includes the kinetic terms for J/'tj) and ip' states and the vertices describing the electromag¬ 
netic vertices XcJ J/P’l and Xcj V'WQ As before we assume the rest frame for the initial state XcJ- The 
kinetic Lagrangian reads 

/:kin(^) = ^2A4 {i{wd) - AM} 4 -)(a:) + ^ 2 M^ {i{u;d) - A'M} V-f )(cc), ( 86 ) 

with the residual masses A A/ = {M‘^ — M^)/2M^ and A'M = (M^ —M^,)/2M^. The fields and 

describes the residual motion of the heavy J/ip and ip' particles and satisfy ui^'tp^^\x) = uj>^ip'^‘^\x) = 0. 
The leading-order in l/m effective Lagrangian describing the radiative decays XcJ J/tp + x ^nd 
-^XcJ + 1 reads 


^SP = Tr k 470^1 + Meg/; Tr k 47o4 


+ h.c. 


with 




?7c75 


7 = +d) {<'“7" - )(i -y^), 


and 


4 = ^(1 + </) |xr7a + Xo + ^(1 


x(l-4- 


(87) 

( 88 ) 

(89) 

(90) 


The currents Js,Js nnd Jp describe particles from S- and P-wave multiplets, respectively. In Eq.(87) 
we introduced the dimensionless couplings and /(,. The fields XJ describe charmonium states XcJ- 


Computing the trace in Eq.(87) one finds 


nem 

^SP 




+eeg 4 xrV'k^^^w'^ -7 xipu:pF^ 


(91) 

(92) 


where we show only the relevant terms. 

^We are grateful to Maxim Polyakov for discussion of the contribution with the virtual state ■>p'. 
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Figure 5: The diagrams which describe the matrix element (96) in HH^PT. The crossed box denotes the 
operator (95), dashes lines describe the Wilson lines and Fr, black circle corresponds to the interaction 
vertices generated by £|p (87). 


Our calculations involve operators 0{^Pj) and ©(^S”!) which have also to be matched onto physical 
quarkonium fields. The spin symmetry in the heavy quark limit yields 




(93) 

(94) 


where af3 are spinor indices. Taking the matrix element and computing the traces one can see that 
Eqs.(93) and (94) reproduce correctly the matrix elements -([^. Using these results one finds 




(95) 


Hence 


(0| \xcj) = (0| T { [{OCS,)) V'(-)"(0) + (O'CS,)) V^'(-)"(0)) Y^Yn, £1’?} \xcj) ■ (96) 

Computing the T-product in Eq.(96) gives the diagrams in F ig{^ These diagrams are UV-divergent and 
we use in our calculation the dimensional regularization and MS subtraction scheme. The results read 


Ixci) = -ie±[aa]e^ i%Q 


(0| e>^('^S'i) \xc 2 ) = i-eq 


where *£_L[tTa] = z£‘^“^^n^np/2 and 
AM 




In 2 — 1 — In 


AM 


tn) + /; (O'^S,)) 


A'M 

M^ 


In 2 — 1 — In 


£x 


-A'M 


(97) 

(98) 

■ (99) 


From these results one sees that the spin symmetry of NRQCD relates the soft photon matrix elements for 
J = 1 and J = 2 which are defined by the same nonperturbative couplings f~^ (Oi^Si)) and /' (©'(^^i)). 
The imaginary part in Eq.(99) corresponds to the photon-quarkonium (J/'0) cut in the diagrams in 
Figj^ Hence we conclude that the two-photon cut appears only in the hard photon contribution. The 
contribution with xj}' has no physical cut and in the diagram this is provided by the negative value of 
A'M. 

The UV-poles in the HHyPT diagrams appear due to the mixing of the operators O'^i^Si) and 
Xcj- Therefore this UV-pole can be absorbed into renormalization of the chiral constant in front of the 
operators Xcj0 The expression for the total amplitude now reads 


- Q! 1 


Al = iUnT2Vn } ■ 


( 100 ) 

( 101 ) 


^Or equivalently one can say that this pole renormalizes the contact vertex describing the Xcj decay 
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We set the value of the chiral scale defining the chiral couplings Cj{no) as product of the 

two-photon hard coefficient functions and constant ^O(^Po)) 

CJ{^,o) = C^4H^io) (Oi^Po))- (102) 

This defines the expressions for the amplitudes which will be used for our numerical estimates. 

4.2 Numerical estimates 


In order to perform numerical estimates we need the values of the nonperturbative parameters 
^O(^Po)) and Two of them are related to the values of quarkonium wave function at the origin, see 
Eqs.([2T|) and (23). Their absolute values have been estimated in Ref. 21 using different models for the 
potential. In our numerical calculations we use the values obtained for Buchmuller-Tye potential |23| 


|P^i(0)|2 ~ 0.075GeV^ 

|Pio(0)|^ ~ 0.81GeV^ |P2o(0)|^ ci 0.530GeV^ 

We also assume that they correspond to positive values: 

.Rio(O) > 0, P 2 o( 0 ) > 0, P 2 i(d) > 0. 


(103) 

(104) 

(105) 


The absolute values of the electromagnetic couplings and /' can be estimated from the decays 
Xcj —> and —)■ XcJl- Using for widths r[xci] = 0.84 x 10“^GeV, r[xc 2 ] = 1-93 x 10“^GeV and 

branching fractions Br[xci —> = 0.340 , Br[xc 2 —> = 0.192 from we obtain 


1 / 7 ! = 


l^[Xcj]Br[xcj -)• J/V’7] 




fc3/M,2 


XcJ 


5.87 (x,o) 
6.05 (x,i) 
6.03 (xc2) 


~ 6 . 0 , 


(106) 


where kg = — M|)/2M^^^ is the photon energy. Similarly, using widths r['0'] = 0.299 x 10 ^GeV 

and branching fractions —>• Xal] = 0.096 , Pr['0' —>■ Xc 27 ] = 0.091 from we obtain 


i/;i = 


6- 5 (Xco7) 

7- 0 (Xci7) 

8- 1 (Xc27) 


~ 7.2, 


(107) 


We also need to know the sign of this coupling which can only be defined by a specific nonperturbative 
calculation. It turns out that this coupling can be represented as an overlap integral of the radial wave 
functions. Comparing our results for the decay amplitudes Xcj J/4’1 with the ones computed in 
Ref. 


24 we find 


1 r°° 

^ = drr3p2i(r)Pio(r), 


(108) 


where factors x'2M appear due to relativistic normalizations of the hadronic states. The analogous 
expression also holds for the coupling /(. The overlap integral has been computed in the framework of 


potential models, see e.g. Refs. 


25 


26 


Its value is found to be positive for P and negative for /'. 


Therefore we assume in the followmg^at /^ > 0 and /' < 0. 
The expressions for the decay width read 


r[xcj 


e~*'e 


^ IC'^V(a^o) {0{^Po)) + C^^eQh{fio)lV2 

(Oi^Po)) + P 


1 

4077 


(109) 


where C-y is given by Eq.([6T|). We use rric = 1.5 GeV for the mass of the charm quark and compute /i(/ro) 
by substituting in the expression (99). 


Our numerical results are presented in Table for different values of iiq. The subscripts s and h 
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Table 1: Numerical results for the decay widths for different values of the factorization scale /xq. 


/xo,MeV 

r[Xci e"e+], eV 

r[Xc 2 ^ e-e+], eV 

300 

0.060s + 0 . 009 , 1 s + 0.023,* = 0.091 

0.036s + 0.020,*s -k 0.016,* = 0.072 

400 

0.063s + 0.013,*s + 0.011,* = 0.087 

0.038s + 0.017,*s -k 0.013,* = 0.068 

500 

0.066s + 0.011,*s -f 0.004,* = 0.082 

0.040s -k 0.015,*s -k 0.010,* = 0.065 


denote contributions from the soft and hard photon terms and hs corresponds to the interference of these 
contributions. In all cases the largest numerical contribution is provided by the ultrasoft matrix element. 
This contribution is relatively large and it weakly depends on the factorization scale /Xq. Our estimates for 
r[Xci e“e^] is approximately factor 5 smaller then the estimate in Ref. and in a good agreement 
with the estimate in Ref. [^. For T[xc 2 e“e“'"] our result is five times larger than one obtained in 
Ref. (^. 

From Table [T] one can observe that the interference of the hard and ultrasoft contributions is numeri¬ 
cally large for Xc 2 width and relatively small for Xci ■ This can be explained as following. The imaginary 
part of /i(^o) is numerically much larger than the real one, see Eq.(99). Further, the hard coefficient func¬ 
tion is real and therefore corresponding interference in the width depends only from the real part of 
/i(/xo). The imaginary part of is not zero and therefore in this case the interference depends on the 
large imaginary part h(/xo) and turns out numerically large. This observation allows one to conclude that 
the decay width Xc 2 is quite sensitive to the relative sign of parameters i?io and i? 2 i- o^r estimate we 
used that these parameters has the same sign, see Eq.( 105). However if they have opposite sign then the 


e e+] by factor 


interference contribution is negative and this reduces the numerical value of the r[xc 2 

2 . 

In Ref. it was shown that unitarity and analyticity allows one to constrain the minimal values of 
decay widths 


r[Xci -t e e+] > -t e e+]r[xci -t 'iJ/fp] « 0.046 eV, 

Z Kq 


( 110 ) 


la 


9a2 


r[Xc2 ^ e e+]> \ Y^bCc2 77 ] + \l 7 ^^[Xc 2 ^ e e+] I « 0.037 eV. (Ill) 


20fc, 


In the presented formalism these constrains are always satisfied because the soft contribution has a 


cut which yields the imaginary part required for the saturation of the bounds in Eqs.(llO) and (1111. 


Therefore all our estimates shown in Table I are in agreement with these inequalities. As one can see 
from Table I the hard two-photon contribution is always smaller than the limiting value in both cases. 
The same observation was also made in Ref. . This clearly indicates that the soft photon configuration 
provides a critically important contribution to these decay amplitudes. 

The derived approach can also be used for a description of leptonic decays of bottonium states XbJ- 
These particles have almost the same branching fractions for Xbj T(15')7 decay, see e.g. 


22 , but at 


present the widths of these states are not yet measured. Therefore, we cannot extract the decay coupling 
using experimental data. Instead, we use the estimates for the corresponding widths obtained in the 

The corresponding values can be found in Ref. 26 and read 


f(b) 

J~t 


model with a Cornell potential 

r[xbi] = 27.8keV, r[xh 2 ] = 31.6keV. 
This gives for the dimensionless coupling in the HHyPT Lagrangian 

4 '^ ~ 9 . 4 . 


On the other hand, the width of T(2S') and branching fractions T(2S') —)■ XbJl are known 22 


Fr[T(25)] = 32keV, Hr[T(2^)] ^ Xbix] = 0.06, Hr[T(2^)] ^ Xb 2 X] = 0.07. 

Using this values we obtain 

~ -16. 


( 112 ) 

(113) 

(114) 

(115) 
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The sign of the couplings f!^'^ and in Eqs.(113) and (115) is again defined with the help of of the 


overlap representation as in Eq.(108) and corresponding estimates given in Ref. 26 . The corresponding 
radial wave functions at the origin read [21| 


|R^i(0)p ~ 2.067GeV®, (116) 

|Rio(0)|^ ~ 14.05GeV^ |i?2o(0)|^ 5.7GeV^ (117) 

With these values and taking /tq = 400MeV we obtain 

r[xw ^ e-e+] = (2.0« + + 1.1,,) x lO'^ - 4.0 x IQ-^eV, (118) 

r[xh2 ^ e-e+] = (1.2, + 0.15/,, + 0.4,,) x lO'^ ~ 1.7 x lO'^eV. (119) 


We observe that in this case the contribution of the ultrasoft configuration also remains larger than the 
hard one. 


5 Conclusions 

The decay width r[xcj e“e^] was computed using a factorization NRQCD approach. The dominant 
partonic subprocess was described by the annihilation of the heavy quark-antiquark pair into two photons: 
cc —> 7*7* —> e“'"e“ . The corresponding contribution is given by the one-loop diagram with two photons in 
the intermediate state. The dominant regions in the loop integral are associated with two configurations: 
hard photons and one ultrasoft and hard photons. The soft part of the contribution with ultrasoft photon 
overlaps with the higher Fock state \QQ^) of the heavy meson, while the hard contribution overlaps with 
the leading two quark state. We have demonstrated that these contributions can be factorized and 
described by two different operators in NRQGD effective theory. The ultrasoft photon contribution is 
estimated using framework of the heavy hadron chiral perturbation theory. This allows us to obtain 
numerical estimates using a minimal set of the known nonperturbative parameters. Our estimates for 
charmonia Xci and Xc2 show that the ultrasoft photon configurations provide the numerically dominant 
contribution. This explains why the obtained numerical results for r[xci e“e“'"] are in good agreement 
with the estimates obtained in Ref. where only the usoft contribution was considered. We also expect 
that the developed formalism can be helpful to perform a more systemic analysis of decays if charmonium- 
like state such as X{3872) —>■ e“e+. 
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